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Abstract. With the present and soon-to-be breakthrough of optical interferometry, countless objects shall be within 
reach of interferometers; yet, most of them are expected to remain only marginally resolved with hectometric 
baselines. 

In this paper, we tackle the problem of deriving the properties of a marginally resolved object from its optical 
visibilities. We show that they depend on the moments of flux distribution of the object: centre, mean angular 
size, asymmetry, and curtosis. We also point out that the visibility amplitude is a second-order phenomenon, 
whereas the phase is a combination of a first-order term, giving the location of the photocentre, and a third-order 
term, more difficult to detect than the visibility amplitude, giving an asymmetry coefficient of the object. We then 
demonstrate that optical visibilities are not a good model constraint while the object stays marginally resolved, 
unless observations are carried out at different wavelengths. Finally, we show an application of this formalism to 
circumstellar discs. 
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1. Introduction 

Optical systems are fundamentally limited in angular res- 
olution by their spatial extent because of diffraction. Very 
soon, the idea of combining light coming from two "dis- 
tant" telescopes was investigated in order to overcome th e 
limitation in size of single pupils. iMichelsonI lll89llll92(t : 
iMichelson k Peasd l|l92l|) derived the angular diameters 
of some solar-system bodies and stars by measuring the 
contrast of the fringes (called visibility amplitude) ob- 
tained when interfering light comes from two apertures: 
this contrast is maximum when these apertures are clos- 
est and decreases with the distance between telescopes 
(called the baseline). The baseline at which the fringes 
disappear holds information on the angular extent of the 
source. After some time, the mid-s eventies saw the come- 
back o f optical in t erfero metry with lHanburv Brown et alJ 
lll974) : iLabevrid (Il975l) . but it long stayed confined to 
bright and simple objects, mostly stellar diameters and 
multiple systems. It is all the more frustrating as the the- 
ory of interferometry allows image reconstruction and as 
radio arrays achieved this goal within a few decades of ex- 
istence: the atmospheric turbulence and the nature of light 
both lead to complex optical designs and have slowed the 
development of optical interferometry. Moreover, the shift 
of the fringes (called phase) is completely blurred by the 
atmosphere, so techniques to retrieve phase information 



— necessary for imaging capacities — needed additional in- 
vestigation. 

Recently, the Pa lomar Testbed Interferometer (PTI, 
IColavita et alJll999l) and the In f rared Optical Telescope 
Array (IOTA. ICarleton et alJ Il994() allowed us to 

[ robe circurnstellar matter in star-forming regions 
Malbet et al.lll998l: lAkeson et allliooot iMalbet fc Berged 
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l2002bHAkeson et al.i200^ giving some constraints on the 
geometry of these objects. With the advent of the Very 
Large Telescope Inte rferometre (VLT I. lGlindemann et alJ 
2000) and the (KI, IColavital l200l|) we are expecting a 
much higher accuracy with their large pupils (8-10 m), 
and good constraints on objects thanks to the num- 
ber of baselines available and partial phase informa- 
tion; yet they will not allow direct image reconstruc- 
tion very soon because recombination will be first per- 
formed with two or three telesco pes. The Navy Prototyp e 
Optical Interferomet re f NPOI. [ Armstrori g et al.' '199^, 
the CHARA array ijten Brummelaar et al. l200Qi) . and 
the Camb ridge Optical Apert ure Synthesis Interferometer 
(COAST, iHaniff et allbood ^ provide imaging capacities 
with a multi-telescope recombination, but with a lower 
sensitivity that renders faint object science difficult. 

We are clearly entering a phase in which more 
than an apparent diameter is measured but no imag- 
ing is performed; in this context, observers and mod- 
ellers use i nterferometric observations as a constraint 
on models (iM albet fc Berged l2002bl lAkeson et all 120021: 
iLachaume et al...200L e.g.), but their link with the ge- 
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ometry of the object remains unclear; it is still quite 
common to think in terms of diameter. For instance, 
iMonnier fc Millan-Gabetl l)2002|) link the uniform disc 
equivalent diameter derived from the IR interferometric 
observations of young stellar objects with the physical ra- 
dius of their supposed inner hole. The phase also raises 
problems of geometrical interpretation. Since it is blurred 
by the atmosphere, one uses the closure phase to retrieve 
partial information on it^: the principle is to add the 
phases over a triplet of baselines provided by three tele- 
scopes, which allows one to cancel atmospheric terms. It is 
generally used either in image reconstruction, at NPOI for 
instance, or as a model constraint. Geometrically speak- 
ing, it is a diagnosis of asymmetry, but its accurate mean- 
ing is seldom made clear enough. 

In this paper, we connect the visibility amplitude and 
phase with the geometry of the object, which allows us to 
retrieve information in a model-independent fashion. In 
Sect. 12 we establish a series development of these quan- 
tities involving the moments of the flux distribution, the 
first ones being the location of the photocentre, the spatial 
extent (diameter), and the asymmetry coefficient (skew- 
ness). It appears as a generalisation of the widespread di- 
ameter measurement. We then draw the consequences of 
the formalism in terms of observation and modelling. In 
Sect. 01 we apply this development to circumstellar discs 
with two examples: the case of an object characterised by 
more than one diameter (star, thermal flux and scattered 
light) and the measurement of the radial temperature law 
in these discs. 

2. Visibility and object geometry 

The Zernicke-van Cittert theorem links the complex visi- 
bility V to the normalised flux distribution / of the object: 



V{u) =11 /(a) exp (— 27rm-Q:) d a. 



where u is linked to the projected baseline B and the 
wavelength A by m = \~^B^ and a the angular location 
on the sky. Optical interferometry usually deals with the 
square amplitude \V\'^ and the phase </>, given by 



tan0 = (ImF)/(ReV^). 



(2a) 
(2b) 



In interferometry two extreme cases are usually dealt with: 
on the one hand, an object is fully resolved if its angular 
size is of the order of B/X as it would be with a single 
dish telescope of size B. In that case, the visibility is arbi- 
trary and highly depends on the shape of the target. On 
the other hand, a point-like source is not resolved and its 
visibility is y = 1. Between these two cases, the object is 
said to be marginally resolved and V is close to unity; its 

^ Other techniques are also the use of a close object as a 
phase reference or the differential phase with spectroscopy, but 
we shall not discuss them in this paper. 



size is a fraction of B/X. In such a case, most of the flux 
is located in a zone where |m-q;| ^ 1, so we carry out a 
series development to derive the real and imaginary parts 
of V: 

Re V{u) = jj I{a) (l - 2iT'^{u-af) d^a, (3a) 

IuyV{u) = jj I{a) ^27r(M-a) - ^TT^{u-af^ d^a. (3b) 

These expressions drive us to define the n-th moment of 
the flux distribution as a symmetric tensor unambiguously 
defined by 



Mn-Ui ■■ - Ur. 



/(a) (a -Ml) ■ ■ ■ [oL-Un) d^a 



(4) 



for any set of vectors iti, m„. (The expressions for 
the first two moments are given in a Cartesian frame in 
Appendix El- With such a definition, the real and imagi- 
nary parts of V become 

ReV{u) = 1-2'k'^M2-u-u, (5a) 
4 o 

ImV {u) = 2t:Mi-u~ -■n-^Mj.-u-u-u, (5b) 
o 

which, after a few calculations, give 

\V\^ = l-ATi'^[M2-u-u-{Mi-uf], (6a) 

</) = -27r(Mi-M) + ^TT^ [Mg-M-lfM 

o (ob) 
-■i{Mi-u){M2-u-u) + 2{Mi-uf] . 

These formulae are quite annoying because the visibil- 
ity amplitude apparently depends on the location of the 
photocentre of the object, given by the first moment Mi 
(see Appendix , and therefore on the pointing accu- 
racy of the instrument. In order to cancel this apparent 
dependence, we define the moments in respect to the pho- 
tocentre of angular location olo by 



(1) M'-Mi---M„ 



(7) 



I{a.){OL — CKq) -Ml 

• • • (a — do) -Mn d^a 
and perform a frame change on Eqs. (|6al6b|) to obtain 
\V\'^ = 1-Ati^{M'2-u-u) + 0{u'^), (8a) 
<P= -2ti {Mi-u) + \tt^ {M'y,-u-u-u) + 0{v!^). (8b) 

o 

In some applications, a higher-order development is 
needed; it is given in Appendix IbI 

Since the phase is not directly measured because of at- 
mospheric turbulence, the closure phase is used instead. 
With three telescopes labelled 1, 2, and 3 simultaneously 
providing the baselines Ui2, M23, and U31, the closure 
phase reads 

= (?!)(mi2) -t- 0(m23) + 0(m3i)- (9) 
Applying Eq. I8bl we derive a concise expression: 

- 4 



-y (M^-Mi2-W23-1t3l) + C(W^)- 



(10) 
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Vij\^ = 1 - ATT^ufj / / I{a,^)p^ dad/3 



D2 

(j) = ^TT^ U12U21U31 jj I{a,P)P^ dadp 

«3 SD3 

Fig. 2. Link between the flux distribution I{a, (3), the visibility amplitude and the closure phase (j) for a marginally 
resolved object. The three ahgned telescopes, numbered 1 to 3, and the basehnes are projected onto the sky. Third-order 
terms and lower have been kept. 



2.1. Resolving an object: size, asymmetry and curtosis along i; they respectively are 



We consider three aligned telescopes labelled from 1 to 
3, providing the baselines M12, M23, and 1431 in a direc- 
tion given by a normal vector i, as represented in Fig. [5] 
The main characteristics of an object we shall consider are 
its mean diameter, its asymmetry and its curtosis defined 



K 



M'^iii 



(11) 
(12) 

(13) 

(M^-i-i) 

They are of common use in statistics. Figure ^ displays 
their link with the shape of a one dimensional distribu- 



(M^ 
Mi 



^3/2' 
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D (mas) D (mas) precision (%) 

Fig. 3. Variation of the square visibility amplitude and closure phase with the angular size of a marginally resolved 
object for a 100 m baseline in K, and the detectability of these quantities. Left panel: square visibility amplitude vs. 
angular size; middle panel: closure phase vs. angular size; right panel: minimum object size needed to detect either the 
spatial extent or the asymmetry vs. the instrumental precision on \V\'^. On the first left two view graphs the detection 
levels for 0.5% and 5% accuracy on \V\'^ measurements are displayed. We assumed a fairly asymmetric object with 
S = 0.5, as weU as = Ay. 



tion: the diameter D is the square root of the variance, 
giving the mean spatial extent; the skewness S increases 
with the asymmetry of the distribution and is zero for a 
symmetrical one; K indicates whether the flux is concen- 
trated in the peak of the distribution or in its wings and 
is zero for a normal distribution. With these notations, 

\V{u)\^ ^ 1 - An^Duf + \ix^ {{K + ^){Duf \ , (14) 

i>^^iT^S{Du)\ (15) 

where the mean baseline u is given by the geometrical 
mean of the three baselines: 

U = \/(Mi2-i)(M23"i)(M31-i)- (16) 

A summary of these results, in a less formal way, is given 
in Fig. 121 

The main implication of these results is that the visibil- 
ity amplitude drop is a second-order phenomenon (YP"u}) 
while the closure phase is a third-order one [SD^u^^. As a 
consequence, closure phase is much harder to detect than 
the visibility amplitude in a marginally resolved object. 
FigureOldisplays the profile of the visibility amplitude and 
closure phase as a function of the baseline for a marginally 
resolved object with a high asymmetry S = 0.5, as well as 
the minima of detection for these quantities as a function 
of the measurement accuracy. It appears that the asymme- 
try is detected for angular sizes 3 to 6 times larger than 
the spatial extent or — which is equivalent — for 3 to 6 
times larger baselines. 

2.2. Validity of the approximation 

The above development presents two limitations: on the 
one hand, it assumes that all moments are defined and, 
on the other hand, the first few terms of the series are no 



longer a good approximation when the object is resolved 
enough. 

In the case of a power law distribution I (a) oc a"', of- 
ten encountered when scattered light dominates, the high- 
order moments are not defined. Therefore the above devel- 
opment is no longer valid. For instance, *Lacha ume et all 
(2001) have shown that a disc with scattering presents 
a quick drop of the visibility amplitude near the origin 
u = 0, that definitely does not present the smooth pro- 
file \V\^ = 1 - 4{TrDu)^. In Sect. O we shaU see how 
to treat scattering at a large scale, while using the above 
formalism for other flux contributions. 

Another important point is the range of validity. The 
left panel of Fig. ^ compares the exact visibility of differ- 
ent types of objects with the second-order estimate (first 
terms of Eq.^J) as a function of the baseline: as expected, 
the approximation is correct for under resolved objects 
with y « 1 but gets poorer with larger baselines. The 
validity of the approximation depends on whether a dif- 
ference can be made between the estimate and the exact 
value, in other terms whether the accuracy of the mea- 
surement is better than the precision of the estimate. The 
right panel of the same figure displays the visibility at 
which the instrument accuracy allows us to measure the 
deviation from the estimate; it appears not to be much 
dependent on the geometry of the object. With a typical 
2% accuracy on \V\'^, the estimate is valid for \V\'^ > 0.8. 

Figure [SI is similar to Fig. ^ but for the closure phase. 
For a typical binary with A. 2, it displays the exact closure 
phase and the third order estimate given in Eq. (|15ll as 
a function of the baseline. The right panel indicates the 
visibility at which a differentiation between them can be 
made as a function of the instrumental accuracy. With the 
typical 2% precision on \V\'^, the estimate of the phase is 
valid for \V\^ > 0.45. 

When the accuracy allows us to see the deviation from 
the third-order estimates, some following orders can be 
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binary 

ring 

gaussian disc 



baseline (arbitrary unit) instrumental accuracy (%) 

Fig. 4. Comparison between the exact visibility amplitude of an object and its second-order estimate. Different geome- 
tries have been assumed to show that little dependence is found: a symmetrical binary (dashes), a ring (dash-dot-dot), 
and a Gaussian disc (dots). The left panel displays the second-order estimate (solid line) and the exact visibilities 
of the different objects vs. the baseline. The right panel displays the visibility amplitude for which an observational 
difference can be made between the estimate and the exact value vs. the instrumental precision; it is the limit of 
validity for this estimate. 




baseline (arbitrary unit) instrumental accuracy (%) 

Fig. 5. Comparison between the exact closure phase of an asymmetrical binary (Amag = 1.2) and its third-order 
estimate. Left panel: closure phase vs. mean baseline. Right panel: minimum visibility for which no observational 
difference can be made between the estimate and the exact phase vs. instrument accuracy. 



measured and the subsequent moments of the distribution 
can be accessed. This can be a means to retrieve model- 
independent spatial information. 

2.3. Model-fitting of observations 

Theoretically speaking, the knowledge of visibility and clo- 
sure phase at all baselines smaller than B would allow 
image reconstruction with an infinite resolution (by us- 
ing the analycity of the visibility); so, marginal resolution 
should not be a problem. This point presents a small but 
irretrievable flaw: it assumes that there is no noise in the 
data. When the finite precision of measurements is taken 
into account, the accuracy on the reconstructed image or 
model is impaired and also, the number of parameters 
constrained in a marginally resolved object. If the visibil- 



ity amplitude alone is available, only the quadratic form 
~ 1 ~ A'K^M2-u-u can be accessed, because the de- 
viation from this law cannot be measured, as shown in 
Sect.O 

We first consider an observation carried out at a sin- 
gle wavelength and use the coordinates (m, v) of u in a 
Cartesian frame. The visibility amplitude and phase then 
are 

V ^l- 4n'^{M'lu^ + 2M'\uv + M'lv^), (17) 

(t)=l - ^TT^{M'lu^ + SM'lu^v + SM'luv^ + M'lv^), 

(18) 

where Af'^ is a component of the rt-th moment Af,'. 
Its expression in the Cartesian frame can be found in 
Appendix El If visibility amplitude alone is considered. 
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Fig. 6. Point-like source model fitting: number of sources constrained by interferometric observations of a marginally 
resolved object vs. number of wavelengths used. Left panel: the fluxes at different wavelengths are not correlated. 
Right panel: the sources are assumed to emit a black-body spectrum. 



the system can be described with three parameters (the 
M'2); a model able to fit any second-order moment will 
fit any data. As a consequence, a marginally resolved ob- 
ject observed in visibility amplitude, whatever the baseline 
coverage is, can be modelled by 

— a ternary system of unresolved stars, 

— a Gaussian elliptic disc, 

— a uniform stellar photosphere and an unresolved star. 

If closure phase information is available, then four ad- 
ditional parameters describe the object (the M'\). A 
marginally resolved object observed in both visibility and 
closure phase can be modelled by a system able to repro- 
duce the second- and third-order moments, that is 

— a system of four unresolved stars, 

— a Gaussian elliptic disc and an unresolved star, 

— a uniform stellar photosphere and two unresolved 
stars. 

It appears that interferometry does not allow us to dis- 
entangle quite different scenarios when the object is 
marginally resolved. This fact is well-known by observers: 
under resolved observations cannot distinguish between a 
uniform and a limb-darkened stellar dis c; neither could 
the first observations of FU Ori by iMalbet et alJ l)l998tl 
exclude either the disc or the binary scenario. As a rem- 
edy, one can use the technique in combination with other 
types of observations, as another model constraint. 

Nevertheless, a multi-wavelength interferometric ap- 
proach can bring more constraints. We consider a A^p 
point-like source model to be fitted to visibility ampli- 
tudes at Nx wavelengths. The fluxes and locations of the 
A^p — 1 first sources constrain the last one, because the flux 
distribution is normalised and centered; therefore there 
are (A^p — 1) locations and Nx{Np — 1) fluxes, that is 
{N\ + 2){Np — 1) free parameters. Observations provide 
3A^A moments of the flux distribution. The characteristic 



number of point-like sources constrained by the measure- 
ments is given by the equality between the number of free 
parameters and that of moments, so that 



JNx_ 



1 



(19) 



The left panel of Fig. El displays the number of sources 
constrained as a function of the number of wavelengths ac- 
cessed, either observing visibility amplitudes only or both 
amplitudes and closure phases. It appears that amplitudes 
and phases can constrain up to 6 point-like sources with 
A^A ^ 10. If we now consider that the sources emit a 
black-body spectrum — or whatever spectrum determined 
by a temperature and a bolometric flux — there are only 
2(A'p — 1) fluxes. The number of point-like sources con- 
strained then becomes 



3^A 



1 



(20) 



The right panel of Fig. shows that the assumption of 
a black-body spectrum allows us to constrain as many 
point-like sources as wanted, provided that the number of 
wavelengths is large enough. For instance, Nx — 10 allows 
us to constrain 9 sources with the visibility amplitudes 
and 14 if phases are also available. 

In conclusion, a multi-wavelength approach allows us 
to fit more parameters and therefore to distinguish be- 
tween different scenarios (disc, multiple system, etc.) 

3. Application to circumstellar discs 

Circumstellar discs are a good target for interferometers 
because they scale from a few tens of AU (where they 
mostly emit thermal light in the infrared) to a few hun- 
dreds of AU (where they present scattered light in the 
infrared) at a typical distance of 150 pc or more. Two is- 
sues of interest are: observing their thermal light, because 
it comes from the first AUs from the star where planets are 
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Fig. 7. Visibility curve of an accretion disc. Solid line: all 
contributions; dashed line: thermal emission of the disc 
only. 

supposed to form, and deriving their radial temperature 
law, because it appears as a good diagnosis of the phe- 
nomena involved (irradiation, flaring, viscous dissipation, 
etc.). In Sect. I8.ll we show how to take into account both 
the thermal and scattered light, which happen to present 
different interferonietric signatures. In Sect. 13.21 we estab- 
lish a connection between the temperature law and the 
wavelength dependence of the visibility. 

3.1. Stellar, thermal and scattered light 

Describing an accretion disc as a marginally resolved ob- 
ject is inaccurate because thermal light occurs at a large 
scale and accounts for up to 10% of the total flux. In order 
to keep the above formalism, we split the image into three 
components: stellar contribution, thermal emission of the 
disc, and scattered light. To each contribution, one can 
associate a corresponding visibility: 



1 if u = 
otherwise. 



(21) 
(22) 

(23) 



where is the mean diameter of the star, Z?th that of 
the thermal emission of the disc, and 5th the skewness of 
distribution of the thermal emission. 

We assumed that both the star and the thermal emis- 
sion of the disc are marginally resolved, that the scatter- 
ing emission is fully resolved as soon as the baseline is 
non-zero, and that the star is too small and symmetric to 
present a phase. Within the approximation that all com- 
ponents have the same photocentre, we derive the total 
visibility 



F,V, + FthVth + FsVs 



where Ftot is the total flux, F^, that of the star, i^th that of 
the thermal emission, and Fs that of scattered light. Fig.[7| 
presents a schematic view of the visibility curve for an 
accretion disc and compares it to that of the thermal light 
alone. It appears that one point of visibility is not enough 
to derive the diameter of the disc, as authors usually do, 
when either the star or scattered flux are present. 

3.2. Temperature profile 

From visibilit ies at different wavelengths, 

iMalbet fc BerseA ^|20n2a^ showed that the tempera- 
ture profile of an accretion disc can be derived. In the 
context of a massive disc, the flux is dominated by 
thermal light so that 



1-|1/(A)|2 = i2iTD{X)B/Xy 



(25) 



To solve this problem, we need a link between the disc 
extent D{X) and the temperature profile. The disc extent 
is given by the second-order moment M2 determined from 
the radial flux distribution of the disc 

F{r) = {T{r)) . (26) 

We chose a self-similar solution for the sake of simplicity 

r(r)ocr"9. (27) 

For a disc presenting an inclination i, the components of 
the moments then write (see Appendix |0 for a demon- 
stration): 



M'l = M'lcos^i, 
M'l = 0, 

so that 

D{X) oc A'/", 
1 - |T/p oc A2/«-2. 



(28a) 
(28b) 
(28c) 

(29) 
(30) 



Using two close enough wavelengths Ai and A2, we deduce 
an estimate of q: 



1 



l log[(l-|^i|^)/(l- 
2 log[A2/Ai] 



(31) 



(24) 



This result appears as a particular case of the multi- 
wavelength approach (see Sect. I2.3|l . that allows us to fit 
more parameters than o r ie visib ility at a single wavelength 
does. iMalbet & Bergeil l|2002al) obtained a similar result 
by a more empirical argument: they state that the ap- 
parent diameter of a disc is proportional to the radius at 
which B\{T{r)) is maximal. It actually conceals the self- 
similarity. 

Note that the result no longer holds when the disc 
presents an inner hole. Figure|Hldisplays the value gapp de- 
duced from Eq. lpTT|l with the H and K bands for a typical 



0.6 
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Fig. 8. Apparent temperature law, given by the exponent gapp as a function of the disc inner truncation (left) and 
the stellar radius (right). Solid line: active disc with q — 0.75 and an accretion rate of 3 x 10^'^ M /yr; dashed line: 
flared passive disc with q = 0.5 and the same effective temperature at 1 AU. The star is a black-body at temperature 
5000 K. The two visibilities are taken in H and K. 



FU Ori disc. The parameters of the disc model are given 
in Appendix^ When the inner gap becomes larger than 
a few stel lar radii, the error on q, |q a,nr> — q\ can be larger 
than 0.1. iMalbet fc Berefed l(2002af) find q,,nn ^ 0.64 for 
FU O ri; with a typical value i?* = 6 Rq ((Lachaume et al.l 
l200l|) . we can estimate q « 0.5 from the curves. 

The stellar radius has also an influence because of the 
unresolved stellar flux, yet, it remains small for FU Ori 
discs (see Fig. IHJ. In the case of a T Tauri star, the error 
could be much larger, because the contribution of the star 
to the total flux becomes important. 



4. Conclusion 

We have developed a formalism that connects the visibility 
amplitude and phase of a marginally resolved object with 
its geometry, namely the moments of the flux distribution. 
It can prove particularly useful when constraining models 
that present analytical moments and allows us to retrieve 
model- independent spatial information in all cases. It also 
establishes that the closure phase of a marginally resolved 
source is a third-order term and the visibility a second- 
order one; therefore, the phase is much harder to detect 
than the drop in visibility amplitude. 

From the formalism, we were also able to estimate 
the number of parameters relevantly fitted to interfero- 
metric measurements. Unless observations are carried out 
at several wavelengths and the model assumes a black- 
body-like emission, only a few parameters can be fitted 
to marginally resolved objects, whatever the number of 
visibility points taken: three point-like source with visibil- 
ity amplitudes only, and a fourth one if closure phase is 
also measured. This limitation is removed when the ob- 
ject is more resolved, that is, if the baselines are longer 
or if the instrumental accuracy is increased, which allows 
us to measure the deviation of the visibility and phase 
from their low-order estimates. This work can therefore 



be seen as a plea for larger baselines than the CHARA ar- 
ray provides, or high accuracy with IOTA/IONIC or the 
forthcoming VLTI and Keck. 

We then applied this theoretical work to circumstel- 
lar discs, by separating the star, the thermal emission of 
the disc, and the scattered light, the two first ones be- 
ing well described by their moments. It also allows us to 
derive, with some hypotheses, information on the radial 
temperature law in these objects, even if they are under- 
resolved, but requires that measurements should be taken 
at two or more wavelengths. This can be applied to other 
field. For instance, limb-darkened stellar photospheres can 
be probed even with underresolved targets: the equivalent 
diameter is dependent on the wavelength and one could, 
with an appropriate model, measure this darkening. 

With high precision measurements and/or multiple 
wavelengths one can access a large number of moments of 
the flux distribution, which theoretically allows image re- 
construction. This is clearly a path that one should inves- 
tigate in the near future. As a particular case, we believe 
it is possible to retrieve the radial temperature profile of 
supposedly symmetrical objects, as as been initiated with 
FU Ori. The method could also prove useful to constrain 
the location of stellar spots with high accuracy measure- 
ments: the link between the location of these spots and 
the first order moments of the flux is much clearer than 
the information given by image reconstruction techniques. 
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Appendix A: Moments in Cartesian coordinates Appendix C: Moments of a self-similar disc 

We use a Cartesian frame in which a has coordinates We consider a face-on disc with a radial temperature law 
{a, 13) throughout this appendix. T(r) oc r""?, where r is the angular distance from the 

The first moment is a vector centre. The flux distribution then reads 

M^ = {MlMl) (A.l) c 

F{r) ^ ^, (C.l) 

giving the location of the photocentre in respect to the (exp — Ij 

origin of the frame. 

^ JJ I{a,(3)adad(3, (A.2a) 

Ml = [ [ I{a,f3)(3dadf3. (A.2b) 



where Ci and C'2 are constants. The reduced moments are 
\x,y)x^dxdy, (C.2) 
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F 
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F 



The second-order moment is a matrix -^^2—77 F{^iy)xy'^xdy, (C.3) 



2 - [ ^4 ^4 j (A.3) M';^ - II F{x, y)y' dx dy, (C.4) 



yielding information on the spatial extent of the object. 

Its components are ^here F is the total flux given by 



M°=jjl{a,l3)a'dadl3, (A.4a) p = JJ F{r)2Trr dr. (C.5) 

A/1 = ///(a,/3)a/3dad/3, (A.4b) 



M2 = yy /(a,/3)/3 dad/3. (A.4c) 

The n-th order moment is a tensor of order n contain- 
ing 2" elements. Since it is symmetrical it is given by only 
n + 1 terms, that we note Af^' for < A; < n: 

= ///(a,/3)a"-V''dad/3. (A.5) 



For the sake of simplicity, we introduce the integral 
Is^ I F(r)27rr^ dr. (C.6) 



Appendix B: V and (p series development 

We define 



By switching to polar coordinates in Eq. (|C.2|) + Eq. 
and with the help of symmetry, we derive 

M'° - h/{2h), (C.7) 
M'l = 0, (C.8) 
M'l = h/{2h). (C.9) 



Ji — Mn-u, (B.l) We perform the change of variables u ~ C2/(Ar in 

Jn = M^^^ u - u, for n > 2. (B.2) Eq. to find 



n times 

The visibility amplitude and phase are then expressed as: 



Is oc A"/'?-'*. (C.IO) 



V = l-A7r'j2 + ^{J, + 3Ji) 

-^(Ja-10J| + 15J4J2) (3 3) M'l ^ M'l ^X^^. (C.ll) 

+ ^ ( Js + 35 J| - 56 J3 Js + 28 Je J2) Since M'l = 0, 

-f 0(1*1°), M'°B^~^M'^B^ 

^ = 2^ Ji - — J3 + — ( J5 - 10 J3 J2) ^ 
3 15 



87^'' / , ^ where ~ Xu is the projected baseline. Therefore 

Jg + 2520 J3 J2 J4 + 560 J| - 36 J7 J2 



— ( J7 + 219 J3 J| - 2IJ2 J5 - 35 J3^5 ; 

4^9 (B.4) 1 - oc A2/«-2. (C.13) 



2835 
-84. 



84 J3 Je - 126 J5 + 756 Js^l - 7560 J3JI) .^^i"^ .^^o* ^^^f'-O'i' second-order moment along 

the minor axis is shrinked by a factor cos i, where i is the 



inclination; but the above demonstration still holds. 
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Appendix D: Power-law accretion disc model 

The FU Ori disc model has been determined with an ef- 
fective temperature T{r) — Kr^^. The influence of the 
central gap or of the star highly depends on the con- 
st ant K. In the standard vis cous disc model [q — 3/4) 
bv lShakura fc SunvaevI lll973l) . 



K = 



(D.l) 



where stands for the mass of the star and M for the 
accretion rate. If a deviation from q = 3/4 is observed, 
we should carry out a proper modelling of the involved 
phenomena. For the sake of simplicity, we assumed that 
K is always determined by the viscous temperature at 
1 AU, so that 



r(r) 



K 



(lAU)- 



(t 



AU 



We took 

= IMo, 
M = 3 X 10"^ Mq/jt. 



(D.2) 

(D.3a) 
(D.3b) 
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